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Abstract 

We present a systematic formahsm for the derivation of the kernel of the BFKL equation 
from that of the GLAP equation and conversely to any given order, with full inclusion of 
the running of the coupling. The running coupling is treated as an operator, reducing the 
inclusion of running coupling effects and their factorization to a purely algebraic problem. 
We show how the GLAP anomalous dimensions which resum large logs of - can be derived 
from the running-coupling BFKL kernel order by order, thereby obtaining a constructive 
all-order proof of small x factorization. We check this result by explicitly calculating 
the running coupling contributions to GLAP anomalous dimensions up to next-to-next-to 
leading order. We finally derive an explicit expression for BFKL kernels which resum large 
logs of up to next-to-leading order from the corresponding GLAP kernels, thus making 
possible a consistent collinear improvement of the BFKL equation up to the same order. 
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1. BFKL evolution and GLAP evolution 

The all-order summation of large logarithms of the center-of-mass energy in perturbative 
QCD is an intrinsically intricate problem, because it requires the simultaneous resumma- 
tion of two large scales: the large energy, and the scale which renders the strong coupling 
perturbative. This difficulty is apparent when trying to fully include the perturbative 
running of the coupling in the evolution equation which sums energy logs in perturbative 
QCD, the BFKL equation. In fact, doubts have been raised in the past on the consistency 
of a BFKL equation with running coupling. It is also apparent when trying to include a 
high-energy (i.e. small x) resummation in the standard renormalization-group (or GLAP) 
evolution equations of perturbative QCD: indeed, until recently it was not clear whether 
(beyond the leading-log level) a resummation of energy logs would be possible using the 
GLAP equation, without having to resort to the BFKL equation. 

The perturbative expansion of the singlet GLAP splitting functions in powers of 
as{Q'^) at fixed x is unstable at small x: indeed, the recently computed NNLO correc- 
tions show explicitly that the instability is visible already around x ~ 10""^. Analogously, 
the perturbative expansion of the BFKL kernel in powers of as{Q^) at fixed is unstable 
at large Q^. In fact, as well known, this instability in practice makes the perturbative 
expansion of the BFKL kernel very poorly behaved essentially for all Q^. Fortunately, 
it has now been demonstrated explicitly [1-6] that simultaneous resummation of the two 
large logs is feasible, and cures both instabilities. This double resummation can be accom- 
plished in a transparent way by exploiting the duality [7-9] which relates the BFKL and 
GLAP equations, whereby one shows that, at the leading- twist level, both equations admit 
the same solutions and thus describe the same kind of evolution, provided their respective 
kernels are suitably matched. 

Whereas at the fixed-coupling level this duality is straightforward to prove, it becomes 
rather nontrivial at the running coupling level. In fact, even the validity of running- 
coupling factorization at small x to any given logarithmic order (i.e. the existence of a 
universal evolution kernel which sums all pertinent logs for any given boundary condition) 
is nontrivial. Running coupling factorization and duality to any logarithmic order at small 
X were first shown in Ref. [9]. The approach of that reference, however, besides being rather 
formal, was essentially geared to computation of small x corrections to the GLAP equation, 
given the BFKL kernel as input. The inverse problem, namely the computation of large 
corrections to the BFKL kernel given the GLAP kernel, though in principle possible, 
in practice requires an all-order computation. This computation is required in order to 
calculate stable BFKL kernels, and thus essential for any kind of phenomenology based on 
the BFKL equation. Furthermore, in the approach of Ref. [9] even the determination of 
next-to-next-to-leading In ^ running- coupling corrections to the GLAP kernel obtained by 
fixed-coupling duality from BFKL is rather cumbersome. 

In this paper, we will develop a rather general approach to running- coupling duality, 
based on the observation that the standard duality relations between GLAP and BFKL 
kernels, which are algebraic equations between kernels in a suitable double-Mellin trans- 
formed space, are promoted to operator equations when the coupling runs. We will show 
that these equations can be solved order-by-order by purely algebraic methods in terms 
of the commutators of the relevant operators, without having to actually determine their 
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spectra, i.e. without having to solve high-order differential equations. This will enable 
us to determine running-coupling corrections to duality in either direction (i.e. when the 
BFKL kernel is determined from the GLAP one or conversely) and for either kind of logs. 
We will check the method and display its power by reproducing known results in a com- 
pact and efficient way, and then we will use it to obtain some new all-order results. As a 
byproduct, we will obtain a very direct re-derivation of factorization at small x. 

The paper is organized as follows. In Sect. 3 we will define our notation, recall some 
basic results about duality at fixed and running coupling, and introduce the general formal- 
ism for the representation of running coupling evolution equations as operator equations. 
Then, in Sect. 4 we will derive within the operator approach the small x resummation 
of the GLAP kernel from the BFKL kernel which was already discussed in Ref. [9]. We 
will then derive a compact expression, based on the Baker-Campbell-Hausdorff formula, 
which generates the corrections to fixed-coupling duality to any desired order (and with 
the running of the coupling given to any desired order). We will check it by reproduc- 
ing the next-to-next-to-leading running coupling corrections to duality which in Ref. [10] 
could only be determined using computer algebra. Finally, in Sect. 5 we will construct 
the large In resummation of the BFKL kernels from the GLAP kernels, which requires 
the computation of an infinite set of running-coupling terms: first, we will show that at 
the leading-log level there exists an operator ordering such that running-coupling duality 
relations reduced to the naive fixed-coupling ones, then, we will show that beyond the 
leading log level naive duality is modified by commutator terms, which we will determine 
explicitly up to the next-to-leading level. This will provide the first determination of a full 
coUinear resummation of the next-to-leading-order BFKL kernel. 

2. Duality with fixed coupling and running coupling 

2. 1 . General formalism 

We will discuss evolution equations for a parton distribution G = xg{x, Q^). This, accord- 
ing to the definition of the factorization scheme, can be thought of as the gluon density, 
or as an eigenvector of a two-by-two evolution matrix in the singlet sector. The kinematic 
variables x and can be thought of as the standard DIS variables, or more generally 
the perturbative scale and a dimensionless scale ratio such that < a; < 1: for example, 
for the production of a final state X in hadronic collisions with center-of-mass energy s, 
= (e.g. the Higgs mass), and x = Q'^/s. We will not consider the dependence on 
other kinematic variables, such as transverse momentum and rapidity, i.e. we will consider 
standard parton distributions, for which ordinary coUinear factorization applies. 

In practice, it is convenient to express the parton density as a function of the logs of 
the relevant kinematic variables: 



We define the Mellin transform with respect to either (or both) of the kinematic variables: 



G^GiCt); 



^ = log-, t = log— . 

X jJL'^ 



(2.1) 




(2.2) 
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G{i,M)^ / dte-'^'Gii.t). (2.3) 

^ — oo 

Note that, by slight abuse of notation, we denote with the same symbol the parton distri- 
butions G{N^ t), G{^, M), and G{N, M), although they are of course different functions of 
the respective arguments. In this paper we will mostly deal with the double-transformed 
distribution G{N,M). 

The GLAP equation takes the form 

^G(iV,t)=7(a.(t),iV)G(iV,t), (2.4) 

where a sit) is the running coupling, while the BFKL equation at the fixed coupling level 
takes the form 

^G(^, M) = x(«., M)G{i, M). (2.5) 

If we expand the GLAP anomalous dimension 7 and the BFKL kernel x powers of 
keeping the respective kinematic variable [N or M) fixed 

7(«„iV) =a,7o(A^) + al-fi{N) + ... (2.6) 
x(«.,M)=a,Xo(M) + alxi{M) + (2.7) 

the GLAP equation (2.4) sums N'^^-'^LO InQ^ terms if contributions up to 7^ are included, 
and the BFKL equation (2.5) sums N'^+^LO In ^ terms if contributions up to Xfc ^^re 
included in the respective kernels. 

The fixed coupling duality relation states that if kernels x and 7 related by 

iV = x(«.,7(«.,Ar)), (2.8) 
M = 7(«.,X(«.,M), (2.9) 

then the BFKL and GLAP equation admit the same solution. This relation is straight- 
foward to derive [1,8] from the observation that the leading twist behaviour of G{N, M) is 
determined by its pole in the M, plane, and that eqs. (2. 8), (2. 9) express the position of 
this pole. 

2.2. Running coupling duality 

At the running coupling level, the BFKL kernel becomes an operator in M space. Indeed, 
the (x, Q^) BFKL equation, obtained by inverse Mellin transform from eq. (2.5) has the 
form 

,2 °o 



^Mi, Q') = r <K^{QVk')G{^. k% (2.10) 

n=0 



where 



dQ^ /Q2^ -M 



Xn(M)=/ -^[j^j (2-11) 
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Therefore, if in eq. (2.11) we let the couphng run with Q^, ag — > q:s(Q^), upon MeUin 
transformation we get 

J oo 

-G{i, M)^Y. &lXniM)Gi^, M), (2.12) 

where ocs is the operator obtained by the replacement t — > — gfj in the expression for 
as{Q^). For example, at the leading log level, where as{Q^) = 0:^/(1 + Poagt) we get 

= . ^ (2-13) 

where ag = q;s(//^). Note that this implies 

[a;\M] = -Po. (2.14) 

If we start again with eq. (2.10), but we let as — > q;s(/c^), upon inverse Mellin transform 
we get instead 

-G(C, M) = 5] Xn(M)a:G(^, M). (2.15) 

Of course, in general as and Xn(-^) do not commute, and therefore the explicit form 
of the functions x„(M) will depend on the argument of the coupling. However, given 
X„(M) which correspond to one ordering, those which correspond to a different ordering 
can be determined by simply computing the relevant commutators. Of course, this is 
just the M-space form of the trivial observation that given the expansion of the kernel 
K{as{Q'^)-iQ'^ /k'^) in powers of as(Q^) we may determine its expansion in terms of agik"^) 
by expressing as{Q'^) in terms of as{k^) through the renormalization-group equation. 
Therefore, at the runnning coupling level the BFKL kernel xic>^s,M) is an operator, and 
its explicit form depends on the operator ordering. Note that if the evolution equations are 
written in double-Mellin (M, A^) space, the GLAP anomalous dimension also becomes the 
operator ^{as,N) obtained by replacing as{t) — > ctg in eq. (2.6), but since ag commutes 
with a function of N the explicit form of 7(0:5, N) does not depend on operator ordering. 

We will assume henceforth for definiteness that as runs with Q^. In such case, the 
running coupling GLAP and BFKL equations in (M, N) space are 

MG{N, M) = 7(a„ A^)G'(A^, M) + Go{N) (2.16) 
NG{N, M) = x(a5, M)G(iV, M) + Fo(M), (2.16) 

where 7(0:3, A^) and x(q;s,-/W") are respectively defined by the expansions eq. (2. 6), (2. 7) 
with Og cxs (and all powers of o^ to the left), and Gq and Fq are boundary conditions. 
Running coupling duality is the statement that given x there exists a 7 such that the 
solutions to eqs. (2. 16), (2. 17), when transformed back to (Q^, x) space, coincide at leading 
twist, and conversely given 7 there exists a x with this property. This result is highly 
nontrivial in that it relates solutions to high-order differential equations. 
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With the leading-order form of ctg eq. (2.13), we can formally expand the running- 
coupling GLAP and BFKL equations in powers of (3q. This generates a series of running- 
coupling corrections to fixed coupling duality, which we can compute order by order by 
solving both equations iteratively as a series in powers of (3o, and then matching the 
solutions. The fact that this matching is possible is nontrivial, however, because if we 
start, say, from BFKL and determine the dual GLAP kernel, the dependence of the 
solutions must match exactly (up to higher-twist corrections), i.e. to all orders in an 
expansion in powers of t up to the given order in /Sq. So, for instance, at 0{(3q) the 
solution contains terms up to 0{t^) in an expansion in powers of t, which must all be 
matched in terms of a single 0{/3q) correction to the duality relation between kernels. 

In Refs. [8,10] it was explicitly checked that this matching is possible up to order 
/9o, and the corresponding running duality corrections were explicitly given for the de- 
termination of the GLAP anomalous dimension from the BFKL kernel. The nontrivial 
nature of this matching suggested that running duality should in fact hold to all orders. 
In Ref. [9] it was shown that this is indeed the case to all orders in (3o, if one uses the 
leading-order form eq. (2.13) of the running coupling, and the leading-order expression of 
X eq. (2.7), i.e., if the BFKL kernel is linear in ctg- Indeed, in this case it is possible to solve 
the running-coupling (M, A^)-space BFKL equation exactly as a differential equation. It 
can then be proven that in the large limit only the inhomogeneous solution to this 
differential equation survives, which is linear in the boundary condition. Furthemore, its 
Q^-dependence can be determined by saddle-point and shown to exponentiate in terms of 
an anomalous dimension which is thereby determined order-by-order in by the saddle 
point expansion, and seen to reproduce the perturbative result. 

Even though it is reasonable to conjecture that these results are quite general, so that 
factorization and duality hold in fact to any logarithmic order, no general argument is 
available. Furthermore, the saddle-point method of ref. [9] is only applicable if x is linear 
in cts and the latter has the leading-order form eq. (2.13). The perturbative brute-force 
matching can of course be applied systematically, but it rapidly becomes computationally 
very cumbersome: in fact, even the computation of 0{(3q) corrections to duality in ref. [10] 
could only be accomplished by means of computer algebra. Here we will present an ap- 
proach which is entirely general and computationally much simpler, even though of course 
still amenable to computer implementation. 



3. GLAP leading singularities from BFKL kernels 

As apparent from eq. (2.13), the running of the coupling is a leading InQ^ effect, but it is 
subleading In ^, and in fact whereas at leading InQ^ terms to all orders in as(3o contribute, 
to nexf^-to-leading In ^ only terms up to 0((q!s/3o)") should be included. Hence, the deter- 
mination of running coupling corrections to duality is simpler when computing the leading 
and subleading In - corrections to GLAP splitting functions from BFKL kernels than when 
computing leading and subleading InQ^ corrections to BFKL kernels from GLAP splitting 
functions: in the former case, the number of terms which must be determined to each 
logarithmic order is finite, while it is infinite in the latter case. In this section we consider 
the simpler problem, postponing the more difficult case to the next section. 
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The problem is to start from eq. (2.17), with the BFKL kernel x computed to some 
finite order in the expansion eq. (2.7), and determine whether there exists a function 
7(q;s,A^) such that, if one uses it in the GLAP equation (2.16) the solutions to the two 
equations coincide at leading twist. In the fixed coupling case, the solution is given by 
eqs. (2. 8), (2. 9), which in particular imply that if terms up to order a" are included in the 
expansion (2.7), then 7 has the form 



7(a„ N) = -is{N/as) + as-iss{N/as) + . . . . 



(3.1) 



In the running coupling case, the problem was solved pcrturbatively in Ref. [8]. Con- 
sider for definiteness the case in which we only include the leading-order contribution xo 
to the expansion of the BFKL kernel x eq- (2.7), and we let — > eq. (2.13). The per- 
turbative solution then proceeds in four steps: first, ols Eq. (2.13) is expanded as a power 
series in asl^Q-^ and the running-coupling BFKL equation (2.17) is solved perturbatively, 
with a result of the form 



G(7V, M) = 



1 



+ 



^Xo)^/ 



(1 - ^x^Y (1 - ^xo) 



+ 



Go{M) 
N ' 



(3.2) 



Second, the inverse M-Mellin transform G{N, t) of this solution is determined up to the 
desired order. Third, this solution is compared to the expansion of the solution to the {N, t) 
space GLAP equation (2.4) to the same order in ag and linear in t, and the corresponding 7 
is determined, as the sum of 7 determined at the next-lowest order, and a running-coupling 
correction. Finally, it is verified that terms corresponding to all higher order powers of t 
which appear at the given order in a exponentiate, so that the solution to Eq. (2.4) is fully 
reproduced to the given order. 

The exponentiation was verified explicitly up to next-to-next-to leading order in 
refs. [9,10], and the running coupling anomalous dimension was found to be ^{as{t), N), 
where 

7(a„ N) = -is{N/as) + as(3Q^iss{N / a^) + (a,/?o)'A7,,,(Ar/a,) + O{a,(3of, (3.3) 
with ^g{N/as) determined by fixed coupling duality eq. (2.8), i.e. as the solution of 

N = asXoils{N/as)), (3.4) 

and 



A7,,(iV/a,) = -xo 



12 



2X(^ 



M=js{N/a,) 



24x[)' 



(3.5) 
(3.6) 



M=7,(JV/a,) 



where all derivatives are with respect to the argument of xq{M)i which is then evaluated 
at ^s{N/as). 
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We can view these running coupling terms as a modification of x, which gets replaced 
by x^^icts, M) = asXo{M)+a'^^{xi{M) + Axi{M)) + .. . from which 7 eq. (3.3) is obtained 
using naive duality eq. (2.8), with 

Axi=/3o^, (3.7) 

and so on. Of course, if higher order contributions to x eq. (2.7) are included, they also 
contribute to 7^, ^ss, etc., as dictated by fixed-coupling duality. 

This next-to-next-to-leading order exponentiation was shown to follow from an all- 
order argument in ref. [9], were the next-to-lcading correction was also rederived as a 
byproduct through a different, though hardly easier method, based on an asymptotic 
expansion of the exact solution of the running coupling BFKL equation (2.17). Here we 
will show that these results can be obtained in a much simpler and general algebraic setting 
without having to actually solve the running coupling BFKL and GLAP equations, either 
perturbatively or exactly. 

3.1. A model calculation 

Let us start by considering a toy calculation where the GLAP kernel is given by its leading 
singularity, i.e. 7(0:5, A^) = N~^as- In this case, the GLAP eq. (2.16) reduces to 

{M -N-^as)G{N,M) = Go{N), (3.8) 

which, multiplying by NM~^ on the left and noting that [N, a^] = [N, M] = reduces to 

{N - M-'^as)G{N,M) = M-^NGo(N). (3.9) 

Recalling [1,8] that the large behaviour of the solution is determined by the pole in the 
M, N plane, this is seen to reduce to the BFKL-like equation 

NG{N, M) = M-^asG{N, M) + Fo(M), (3.10) 

with a suitably readjusted boundary condition Fq{M) = M~'^6lsGq{M~^6ls). Hence, the 
running dual of 7 = N~^as is 

= asM~^ + PoalM~'^ + 0{al), ^ 

where we have used the commutator 

[as, M-^] = -/3oasM-^as = -/3oM-^alM-\ (3.12) 

which may be readily obtained from eq. (2.14). 

Now we note that the 0{as) term is just the fixed-coupling dual of 7, while the 0{a1) 
term coincides with minus the running coupling correction eq. (3.7): with xo = ^/M, 
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= — This is as it should be, because Eq. (3.5) gives the addition to 7 due 

to the running-coupUng corrections which are necessary for it to be the running couphng 
dual of a given x- Equation (3.7) then expresses this quantity as an effective contribution 
to X- Equation (3.11) instead gives the x which we must start from in order to obtain 
7 = as/N through running-coupling duality. Clearly, this x is obtained by undoing the 
effect of running duality, i.e. by starting off with a x which differs from the naive dual by 
a term which is then removed by the running coupling correction. 

We conclude that in this simple case we have in fact reproduced the running-coupling 
correction to duality. We notice furthermore that the boundary condition did not play any 
role in the argument: in fact, the above argument can be recast as the simple observation 
that the operator equation 

NG{N, M) = M-^asG{N, M), (3.13) 

can be rewritten as 

MG{N, M) = N-^asG{N, M). (3.14) 

The running-coupling duality corrections are then simply due to the commutators which 
are necessary in order to bring the the BFKL kernel in the form of eq. (2.12), which in 
this simple example are straightforward to compute. 

In the next subsection we shall prove that in fact the determination of running coupling 
duality corrections can always be reduced to the algebraic manipulation of an operator 
equation, which we will then perform systematically in the subsequent section. 



3.2. Factorization 

In order to discuss the form of the duality condition, it is more convenient to start from 
the GLAP equation (2.16), which can be formally solved as 

G{N, M) = (M- 7(a„ N))'^ GoiN). (3.15) 

This equation formally embodies all order (GLAP) factorization, in the sense that on the 
right hand side the perturbative evolution factorises from the nonperturbative boundary 
condition Go{N). We can expand this solution about any N as 

G{N, M)=[M- ^{as, N) - 7'(a„ N){N - N) + 0{{N - Nf)] (G'o(iV) + 0{N - iV)), 

(3.16) 

where the prime denotes partial differentiation with respect to N. 

Now, the leading twist behaviour of G{N,t) at large t is given by the pole in the 
M-plane due to the vanishing of the denominator in (3.15): at the pole 

MG{N,M) = -f{as,N)G{N,M). (3.17) 

Clearly the position of the pole will depend on A^. Now assume that there exists an 
operator- valued function x{ois,M) such that at the pole eq. (3.17) may be inverted, i.e. 

NG{N,M) = x{as,M)G{N,M) (3.18) 
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The operator xi^s-, M) is the inverse of the operator 7(ci!s, N), in the sense that 

MG{N, M) = 7(a„ x(«., M))G(iV, M), (3.19) 

and conversely 

NG{N, M) = xias, li&s, N))GiN, M). (3.20) 

Of course the existence of the function x(q;s, M) is nontrivial: we will show how it may be 
constructed order by order in the next section. 

Given x(<^s, -^), we may then choose N in the expansion (3.16) to be the position of 
the perturbative pole in the A^-plane for a given M: using (3.19) the first two terms in the 
square brackets then cancel by construction, and we are left with 

G{N, M) = {N- x(«., M)))-^ (-y X(«., M)))-iGo(x(«., M)) + . . . , (3.21) 

where the dots denote terms which are regular in the neighbourhood of the pole. We may 
then define the function 

Fo(M) = (-7'(a„x(a.,M)))-iGo(x(«.,M)), (3.22) 

which is just some function of M, related in a complicated way to the original function 
G'o(M), but with manifestly no dependence on N . Identifying Fq{M) with a BFKL bound- 
ary condition, eq. (3.21) is the solution to the running coupling BFKL equation (2.17), 
with the kernel x{(^si M) given by eq. (3.19). The solution is again factorised at the pole 
into a perturbative evolution factor and a nonperturbative boundary condition. Thus all 
orders factorization of the running coupling BFKL equation follows directly from all orders 
factorization of GLAP provided only that the kernel defined through eq. (3.19) actually 
exists. 

It is important to understand in which respect eq. (3.19) differs from the naive fixed 
coupling duality relation eq. (2.9), to which it bears formal resemblance. The meaning of 
this equation is that at the pole, the operators and xi^^s^M) have the same action on 
the physical solution G{N,M), even though they are in general different operators. The 
anomalous dimension 7(q;s, N) is then the operator such that if we apply it to a physical 
state at the pole, and then use repeatedly eq. (3.18), we obtain the l.h.s. of eq. (3.19), i.e. 
M. Just as in the fixed coupling case the duality relation (3.19) is invertible: it can be 
used to construct x given 7 or 7 given x- 

To make this more explicit, consider the simple example of sect. 3.1. Starting with 7 = 
7V~^ds, the pole condition (3.17) becomes MG = N~^asG = asN~^G. Then construct x 
by demanding that it satisfies (3.19): 

MG = a.ixi&s, M))-^G. (3.23) 

Note the ordering here: only when the factor of is on the right, and thus able to act 
on G, can we use (3.18). Expanding x at lowest nontrivial order we get 



MG=[as{{e,sXo{M) + alxi) + ---) 

(3.24) 



= [«,(xo(M) + «,xi(M) + ---)"'«j']G' 



[ (Xo(M) + a,Xi(M) + • • •)"' - a./^oll^jC, 
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where in the last hne we used the commutator 



[a-\xo\m = -/3o(Xo ' W)' = ^ox'o/xl (3-25) 

We immediately get that to order o;,, xo(-^) = 1/-^, and substituting this back in 
eq. (3.24) and expanding xi(M) = /3o/M^, in agreement with eq. (3.11). 

Hence, wc sec that the running coupling dual is determined by solving cither of the 
two operator equations (3. 20), (3. 19) from the respective pole conditions (3. 17), (3. 18). In 
the next subsection, we will construct a solution of these equations which, when expanded 
out, generates running coupling corrections to any desired order. 



3.3. A Systematic Algebraic approach 

The construction of the running coupling dual of a given BFKL kernel or GLAP anomalous 
dimension presented in the previous subsection can be formalized in the following way: 
given an operator equation of the form 

pG = qG (3.26) 

(where G denotes the physical solution), and given a function /(g), determine the function 
g such that using eq. (3.26) on physical solutions one gets 

f{q)G = g{p)G. (3.27) 

Note that f ^ g necessarily, since p and q do not commute: for example 

fG = qpG = {pq + [q,p\)G = (p^ + [g,p])G ^ p'^G 

if [QjP] 7^ 0. So the problem is in general nontrivial. 

Equation (3.18) corresponds to the particular choice 



p = aJ^N 



(3.28) 



where for future convenience we have multiplied both sides by as- Equation (3.17) corre- 
sponds instead to the choice 

p= M 

(3.29) 

q = ^{as,N). 

In the latter case, for example, we can choose the function f = Xi such that xili^^s, N)) = 
N: because ojg and commute this is just the naive dual. Equation (3.27) then determines 
the running coupling dual x(M) as the function g{p). 

A general solution to this problem may be obtained by using the Baker-Campbell- 
Hausdorff formula for a pair of noncommuting operators A and B which to cubic order 
is [11] 

e^e^ = exp{^ + B+^[A,B] + ^{[A, [A, B]] + [B, [B, A]]) + ...}. (3.30) 
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Letting A = q and B = p — q we get 

e'^e^"^ = exp{p- ^\p,q] + [q,p]] + ^\p, \p,q]] + ^[g, [q, [q,p]]] + -^[g, [p, [p, g]]] + ....}. 

(3.31) 

Multiplying the right-hand side by the identity e^e ^ = 1 on the left, and using the Baker- 
Campbell-Hausdorff formula again with A = —p and B set equal to the exponent of the 
exponential on the r.h.s., eq. (3.31) becomes 



e^gP-g = ePexp{-|[p,g] + i[g, [g,p]] + [p,g]] + ^[g, [g, [g,_^^ 

I O.O.Z I 

+ [p, \p,q\]] -i\p,\P, \p, q]]] + ■■■}■ 

Now, eq. (3.26) implies that 

eP-'^G = G, (3.33) 

so eq. (3.32) gives us an expression for the action of on physical solutions G. We can 
use this result to construct a general solution to our problem by noting that any function 
/(g) can be constructed by Taylor expansion as 

/(g) = e^A/(A)|^^^. (3.34) 

Rescaling the operators p and g in eq. (3.32) by i.e. letting p — > g — > g^ we get 
finally 



MG = e^A exp{-i[p, g]^ + i[g, [q,p]]^ + [p,^]]^ + ^Q, [q, [qM£ 

ir^ .-m 



(3.35) 

Expanding out the exponential on the r.h.s. of cq. (3.35) leads to an expression in terms of 
f{p) and repeated commutators of p and g. For example, treating the simple commutator 
as first order, the double commutators as second order and so on, expanding up to second 
order we get 

f{q)G = {m - y"{p) [p, g] + y"'{p) [g, [q,p]] 

„ (3.36) 

+ y"'ip) [p, [p, g]] + irm, + .)}g. 

This is the main result of this section. In order to show how it can be used in practice, 
in the next subsection we will use it to reproduce the next-to-next-to-leading order result 
eq. (3.6). 

3.4- Running coupling corrections to NNLO 

In order to compute running- coupling corrections to the duality relation between xo s-nd 
the expansion eq. (3.1) of 7 in powers of ag at fixed ag/N we use eq. (3.36) with the 
identification eq. (3.28) of p and g. Furthermore, we assume that x is given by its leading 
order expression, 

x{as,M) = asXo{M), (3.37) 
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so q = xo{M). The relevant commutators are then 



[p,q] = [a:'N,xo] - -iV/?oXo(M) 

[q, = [xo(M), [a;'N,xo{M)]] = 0, 

where xo and its derivatives are all functions of M. 

We now choose as function f{q) the function 7^, which is related to xo by naive 
fixed-coupling duality eq. (3.4), so 

M = 7,(xo(M)). (3.41) 

Substituting in eq. (3.36) we get 

M = {ls{&;'N)-lyj{a:'N) i-NPox'oiM)) + H'i^s'N) ((iV/3o)^x()'(M)) 

(3.42) 

where the prime denotes difi'erentiation with respect to M of xo and with respect to Q!~-'^A'" 
of 7s. Note that in this equation the first term on the r.h.s. is 0{a'^), the second term 
(coming from the simple commutator in eq. (3.36)) is 0{N) = 0{as) in an expansion in 
powers of as at fixed aJ^N, and the last two terms, coming from the double commutator 
and the commutator square, are 0{N'^) = O(d^). 

The right-hand side of eq. (3.42) does not yet give the sought-for expression for 7, be- 
cause of the residual dependence on M through xo and its derivatives. This dependence can 
be eliminated iteratively by solving the equation to leading order in as, back-substituting 
the result to determine the NLO solution and so on: to LO of course we simply get 

MG = {-fsia-^N) + 0{as)}G, (3.43) 

which replaced in the 0{a1) term gives 

MG = ilsiaj'N) + ^NPol'J{&;'N)x'o bsi&j'N)] + ©(d^)} G. (3.44) 

Beyond 0(6:^) however the back-substitution becomes nontrivial, because an operator 
function 0{M) on a physical solution on which M = 7s (ci^^-^) docs not simply become 
(9(7s(d7^-^)), for precisely the same reason why the functions / and g in eq. (3.27) are 
different. Specifically, to O(a^) the back-substitution of the leading-order result into the 
0{as) term in eq. (3.42) requires the evaluation of 

x'o{M)G = {x'o bs{&;'N)] + ^N(3ol's{&j'N)x'^' bs{&j'N)] + 0{al)}G, (3.45) 

where the result follows using eq. (3.35) again, but with p and q given by eq. (3.29), and 
/ = Xo- Note that it is suflBcient to evaluate eq. (3.45) to 0{as) because the term on the 
left-hand side appears in a contribution to eq. (3.42) which is already 0{as)- Using this 
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result, the back-substitution of the next-to-leading order result eq. (3.44) into the 0{aa) 
term in eq. (3.42) has the form 

(3.46) 

+ x'o'(7s)im7sXo(7s)] + 0(d^)}G', 

where the second term in square brackets comes from eq. (3.45), and the third from the 
expansion of the argument of x' about 7^. All terms on the right hand side of this equation 
now depend on aJ^N through 7^ and its derivatives. 

Collecting everything, the final result for 7 is given by the r.h.s. of eq. (3.42), with 
the 0{as) term evaluated using the next-to-leading back-substitution eq. (3.46), and the 
two 0{a'^) term evaluated with the leading back-substitution eq. (3.43), to give 

MG = {7. + NPohsX'oils) + {NPomilsfxoils)Xoils) + hWsXoils) 

+ H'x'^ils) + bs"ix'oils)r) + Oial)}G. 

In this way, all terms in eq. (3.47) depend on aJ^N through 7^ and its derivatives, while 
the dependence on M has been eliminated. The result can be further rewritten entirely 

in terms of xo and its derivatives, or in terms of 7s and its derivatives, by exploiting the 
fact that 7s and xo and their derivatives are related by the fixed coupling duality relation 
eq. (3.41) and the following relations which may be obtained differentiating it: 

7^ = VXo, 

1': = -x'M' 

2 3.48 

Is = (3Xo - XoXo )/Xq 

1'.:" = -(15xf - W,XW^' + x[,'xrXo)/x[,', 

where xo and its derivatives are evaluated as a function of ^s{,ol~^N) (so xo = oc~^N). 
With these substitutions, eq. (3.47) takes the form 

MG = {ls{&-'N) + apo^^ssi&j'N) + {aPofA^sss{&s'N) + 0{al)}G, (3.49) 

with A7ss and A7sss given by eq. (3.5) and eq. (3.6) respectively, thus reproducing the 
result of ref. [8,10]. It is easy to extend this calculation to higher orders, and also to include 
the contributions due to higher-order terms in the (3 function, which will simply modify 
the form of the basic commutator eq. (2.14). A systematic investigation of such terms is 
presented elsewhere [12,13]. ^ 



^ When going to yet higher orders two further subtleties appear. First, eq. (3.36) cannot be 
used directly, but rather one must use the modified form of it which has the function f{p) all the 
way to the right in all terms. This is due to the fact that it is only when acting on the physical 
solutions on the right that one can let p — q iteratively. The second is that when expanding the 
argument of xo and its derivative about 7s it must be taken into account that 7s is shifted by an 
amount that does not commute with 7s itself, so the Baker-Campbell-Hausdorff formula for the 
function of the sum of noncommuting operators must be used. 
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4. BFKL leading singularities from GLAP anomalous dimensions 

Let us now turn to the inverse problem to that which we discussed in the previous section: 
namely, the construction of a BFKL kernel which is dual order by order to a given GLAP 
anomalous dimension eq. (2.6). Obviously, at the fixed coupling level, duality maps the 
expansion eq (2.6) of 7 in powers of as at fixed N onto the expansion of x in powers of cug 
at fixed ag/M: 

xias,M) = XsiM/as) + (M/«,) + ..., (4.1) 

so that the dual of 7 when terms up to 7^ are included in eq. (2.6) is given by x with terms 
up to Xs" included in eq. (4.1). Now, at the running coupling level, the very definition of 
Xs requires an operator ordering, since ag and do not commute. Moreover, eq. (3.12) 
implies that different operator orderings within Xs differ by terms which are of the same 
order as Xs itself (and similarly to subsequent orders). Hence, as already mentioned in 
section 2, the determination of the dual Xs involves the computation of an infinite set 
of running coupling corrections to duality. Running duality contributions to Xs may be 
computed to any desired order by using the method of the previous section, but with the 
identification eq. (3.29) of p, q. This may be useful for the determination [13,14] of an 
approximate form of Xn beyond the known next-to-leading order xi- However, in this 
section we will show that, by choosing a suitable operator ordering, it is in fact possible to 
determine runnig coupling corrections to all orders in as, and we will construct explicitly 
the running coupling dual to 7 eq. (2.6) up to next-to-leading order. 

4-1- Operator ordering and leading order 

In order to introduce our proof, consider first a simplified case. Namely, assume that we 
start with the model leading-order GLAP anomalous dimension 

7toy = «.(t)(A^-'-l). (4.2) 

This anomalous dimension has the same small-A?" behaviour as the leading-order GLAP 
anomalous dimension (up to an overall numerical factor, which is irrelevant here), and 
it respects momentum conservation in that it vanishes at = 1. In fact, the simple 
form eq. (4.2) turns out to be a surprisingly good approximation to the full leading-order 
GLAP anomalous dimension, even though this is immaterial for our purposes. We wish 
to determine the Xs kernel which is dual to it. Because the anomalous dimension is very 
simple, we can solve all the relevant evolution equations exactly. 

At the fixed coupling level, the dual eq. (2.8) of 7toy is given by 

Xtoy{M/as) = = . ^ I., ■ (4.3) 

as + M 1 + M/as 

We proceed by trial and error: we guess a particular ordering in x, we determine the 
corresponding 7 anomalous dimension explicitly, and we compare the result to eq. (4.2). 
Specifically, we make the ansatz 
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where in the last step we used the commutator eq. (2.14). This can be thought of as 
the result of writing the fixed-coupling dual (4.3) as a function of a~^M, and then let- 
ting as — ^ tts with this particular ordering. Substituting in the running-coupling BFKL 
equation (2.17) we get 

{1-Pq + MaJ^)NG{N, M) = Gq{M) + G{N, M), (4.5) 

where Go{M) = (1 — /?o + MaJ^)G{N, M) is a modified boundary condition. We now 
assume that dig takes the leading-order form eq. (2.13). As proven in Ref. [9], the large 
behaviour of the solution to eq. (4.5) is the same as that of the solution of the associated 
homogenous equation, multiplied by a Q^~independent boundary condition. Hence, we 
determine it setting G = in Eq. (4.5), which then has solution 

G{N,M) = M^o-'(i-^-')-ie^/^o"^ (4.6) 

In order to determine the large behaviour explicitly, we must invert the M-Mellin 
transform: 

^ W 27rz . (4.7) 

= (/3oa.(t))^o-Hi-A^-^)[r - N-') - 1)]-^ 

The anomalous dimension is then 

7(iV,a,(t)) = ^^lnGiN,t) = a,(t) (iV"! - l) , (4.8) 

which coincides with eq. (4.2). Hence, we find the surprizing result that for the particular 
ordering (4.4) the running-coupling dual coincides exactly with the naive dual. 

Let us now prove that this result holds in general at the leading In level, provided 
the operator ordering of eq. (4.4) is adopted. As discussed in sect. 3.2, the dual BFKL 
kernel is constructed by starting with the pole condition eq. (3.17), and determining the 
operator xi^siM) which satisfies eq. (3.20). Now, at leading order the pole condition is 
just 

aJ^MG{N, M) = -fo{N)G{N, M). (4.9) 
Therefore, given the ordinary inverse function Xs 

Xs{lo{N)) = N (4.10) 

then the operator-valued function Xs{ci~^^M), when acting on the solution G, using re- 
peatedly eq. (4.9) is seen to satisfy 

NG{N, M) = Xs{a-^M)G{N, M). (4.11) 

Namely, each time a~^M acts on the physical state G{N,M) it gives 'j{N), so acting 
repeatedly gives powers of the commuting operator 'y{N). Therefore if Xs is chosen to 
satisfy the inverse function condition eq. (4.10) the desired dual BFKL equation eq. (4.11) 
immediately follows, which is what we set out to prove. Note that {aJ^M)~^ = cksM"^, 
so we can equivalently take Xs to be a function of cksM"-*^. 
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4-2. Next-to-leading corrections 

The result of the previous subsection can be derived equivalently by letting 

p = a~^M 

. 1, (4-12) 

so that when 7(0:5, A^) = asjQ{N), then p and q commute, and eq. (3.36) implies that 
f{p)G — f{q)G so taking f = Xs defined by naive duality (4.10) we get 

Xs{p)G{N, M) = Xs{lo{N))G{N, M), (4.13) 

whence the result. However, at next-to- leading order q — 70 (iV) + q;s7i(A^), and \p,q] = 
—Po&s'fiiN), which is of the same order as 71 (AT). Likewise, all higher order commutators 
[p, [p, . . . , IP J q]]] are of the same order a^. Therefore, the iterative commutator approach 
of the previous section does not help and we must proceed to all orders. 

This can be done by noting that the desired next-to-leading order contribution 
Xss{(^^^M) to the BFKL kernel must satisfy [recall eq. (3.19)] 

boiXsiP) + &sXss{p)) + asliiXsiP))] G{N, M) = pG{N, M), (4.14) 

where p is given by eq. (4.12). At leading order, this just gives 

lo{Xs{p))G{N,M)=pG{N,M), (4.15) 

which is equivalent to eq. (4.13). In order to determine the next-to- leading order Xss it 
is necessary to expand out the argument of 70 in eq. (4.14), which is nontrivial because 
[«S7 Xss(p)] 7^ 0. It turns out, however, that this commutator can be determined explicitly, 
and this is sufficient to solve eq. (4.14) by expansion. 
To this end, write Xss as a power series: 

X..(p) = J^/Cip-'. (4.16) 

i 

In practice we will find that only positive powers of i contribute to the sum (since Xss 
vanishes when ag 0), but this is not necessary for our argument. Hence, we wish to 
determine the commutator [p~^, dts]- To this purpose, first we note that the commutator 
eq. (3.12) implies 

[&s,P~^] = -f3oP~^asp-\ (4.17) 

which can be rewritten as 

p-^&s = &sp-^ (1 - f3op-Y' = ^s{p - f3o)-^- (4.18) 
Iterating this equation, it immediately follows that 

p-"«. = «.(p-/5o)-". (4.19) 
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Similarly, 

[a„p]=/?o«. (4.20) 

implies 

r^s = (1 - f^oP'T = ^s{p- PoT, (4.21) 

or, equivalently, eq. (4.19) holds with both positive and negative n. Using this in the 
expansion eq. (4.16) of Xss implies that 

xip)as^asXip-Po). (4.22) 

Of course, this result holds for any function that may be expanded as a series of powers 
of p. 

It is now easy to determine the function on the left-hand side of eq. (4.14), i.e. 

7o(Xs(p) +Q!sXss(p)), by expanding 



7o 



{N)^J29jN'^ (4.23) 



where in practice of course the sum runs over — 1 < J < oo, though again we will not 
actually use this in our argument. We must therefore evaluate 

iXs + asXssT = Xs + (xr'«sXs. + Xs~'^asXssXs + ■■■ + ^sXssXs'^) + 0{al), (4.24) 
where Xs and Xss sixe both functions of p. Using eq. (4.22) we get 

iXs + &sXssT = xliP) + &s ^^f.~^'\~^^l% ss{p) + 0{al), (4.25) 

Xsip- Po) -Xsip) 

where we have made use of the fact that [xsip — Po),Xs{p)] — and we have noticed that 
for a pair of commuting operators A and B, 

An _ r>n 

^ + . . . + = — . 

A- B 

Furthermore, it is easy to see that eq. (4.25) also holds for negative n, by expanding 

iXs + ttsXss)"" = (X7^ + Xl^O^sXssXs^T + ^i^l) 

= xr" + (xr"«.x..x;' + xr^""'^«.x..xr' + • • • + x;'«.x..xr) + o(«'), 

(4.26) 

and then using again eq. (4.22) but with Xs ^ X7^- 

Substituting eq. (4.25) in the expansion eq. (4.23) of 7 gives 

loKXsKP) + oisXssKP)) = loKXsKP)) + otsXssKP) _r\_ — 'F\ (^-^^i 

Xs\P Po) XsyP) 
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Using this in eq. (4.14) with the leading-order eq. (4.15) we get the remarkably simple 
result 

XssiP) = Po'liiXsiP)) {Xsip - Po) - Xsip)) . (4.28) 

This is the running coupling generalization of the well-known fixed coupling result 
Xss{M/as) = - jl[^l[Z^/al)) ^ which it reduces in the limit Po ^ (since x's = l/7o(Xs))- 
The result can be generalized to higher orders and systematized by performing the 
expansion of the higher-order generalizations of operator duality condition eq. (4.14) by 
means of the Baker-Campbell-Hausdorff formula for the Taylor expansion of a function of 
the sum of two non-commuting operators. 



5. Conclusions 

In this paper we have presented an explicit constructive proof of the fact that the leading- 
twist duality between the BFKL and GLAP equations holds to all orders, with the running 
of the coupling included to any desired order. Besides its conceptual interest, this result 
allows one to exploit fully the information contained in fixed-order determinations of the 
BFKL and GLAP kernels by using each to resum the unresummed singularities in the 
other. 

Because of the instability of the fixed order expansion of the BFKL kernel in powers 
of as due to unresummed coUinear singularities, our running duality is especially useful in 
constructing a resummed form of this kernel, using the coUinear resummation implicit in 
the fixed order expansion of the GLAP kernel. The running duality can then be used again 
to derive from the resummed kernel a corresponding anomalous dimension in which small 
X singularities are resummed. Indeed, several of the results contained in this paper were 
used already in the explicit construction of a fully stable resummed perturbative expansion 
of small X anomalous dimensions [2,3,6]. Further investigations of the application of the 
formalism presented here to this and related problems will be presented elsewhere [13]. 
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